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ABSTRACT 



We study the geometry and the physical conditions in the inner (AU-scale) cir- 
cumstellar region around the young Herbig Be star MWC 147 using long-basehne 
spectro-interferometry in the near-infrared (NIR /('-band, VLTI/AMBER ob- 
j^ I servations and PTI archive data) as well as the mid-infrared (MIR A^-band, 

VLTI/MIDI observations). The emission from MWC 147 is clearly resolved and 
has a characteristic physical size of ~ 1.3 AU and ~ 9 AU at 2.2 /im and 11 /xm 
00 ' respectively (Gaussian diameter). The MIR emission reveals asymmetry consis- 

^ ' tent with a disk structure seen under intermediate inclination. The spectrally 

t:j- ■ dispersed AMBER and MIDI interferograms both show a strong increase in the 

characteristic size towards longer wavelengths, much steeper than predicted by 
P ■ analytic disk models assuming power-law radial temperature distributions. 

O ■ We model the interferometric data and the spectral energy distribution of 

MWC 147 with 2-D, frequency-dependent radiation transfer simulations. This 
/\ • analysis shows that models of spherical envelopes or passive irradiated Keplerian 

c3 ■ disks (with vertical or curved puffed- up inner rim) can easily fit the SED, but 

predict much lower visibilities than observed; the angular size predicted by such 
models is 2 to 4 times larger than the size derived from the interferometric data, 
so these models can clearly be ruled out. Models of a Keplerian disk with optically 
thick gas emission from an active gaseous disk (inside the dust sublimation zone), 
however, yield a good fit of the SED and simultaneously reproduce the absolute 



^Based on observations made with ESO telescopes at the La Silla Paranal Observatory under programme 
IDs 074.C-0181, 076.C-0138, and 078.C-0129. In addition, this work is based in part on archival data obtained 
with the Spitzcr Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute 
of Technology under a contract with NASA. 
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level and the spectral dependence of the NIR and MIR visibihties. We conclude 
that the NIR continuum emission from MWC 147 is dominated by accretion 
luminosity emerging from an optically thick inner gaseous disk, while the MIR 
emission also contains contributions from the outer, irradiated dust disk. 

Subject headings: accretion, accretion disks - stars: formation - stars: pre- 
main-sequence - stars: individual: MWC 147 - techniques: interferometric 



1. Introduction 

The spatial structure of the circumstellar material around Herbig Ae/Be (HAeBe) stars, 
i.e. intermediate-mass, pre-main sequence stars, is still a matter of debate. Until recently, the 
spatial scales of the inner circumstellar environment (a few AU, corresponding to ^O'.'l) were 
not accessible to optical and infrared imaging observations, and conclusions drawn on the spa- 
tial distribution of the circumstellar material were, in most cases, entirely based on the model- 
ing of the spectral energy distribution (SED). However, fits to the obs erved SEDs are highly 



ambi guous, and very different models such as spheri cal envelopes (JMiroshnichenko et al. 



19971 ). geometrically thin disks ( Hillenbrand et al 



velope (JNatta fc Kruegellll995 : 



Natta et al 



disks, puffed-up inne r disk rims ( iDuUemond et al 



2001 



19921). disks surrounded by a spherical en- 



Miroshnichenko et al.lll999l ). flared outer 



200ll ). and disk plus inner halo models 



(IVinkovic et al.ll2006l ) have been used to successfully fit the observed SEDs of HAeBes. It was 
also unclear how the accretion of circumstellar material onto the star affects the structure 
and emission of the inner circumstellar environment. 

In the last couple of years, long-baseline interferometry in the near- (NIR) and mid- 
infared (MIR) spectral range provided, for the first time, direct spatial information on 
the inner circumstellar regions of young stars a t scales of < 10 na i lli-arcseconds (mas; 
e^ 



Malbet et al. 



11998 



Millan-Gabet et al.l 



1999 



Akeson et al. 12000: iMillan-Gabet et al. 



2001 



Wilkin fc A 



ceson 



20031: 



Eisner et al.l l2004l : iLeinert et al.l l2004l : iMonnier et al 



2005 



Malbet et al.l 120071 ). The interferometric data were usually interpreted by fitting simple an- 
alytic models for the brightness distribution to the measured visibilities. The characteristic 
sizes of the emitting regions derived in this way were found to be correlated to the stellar 
luminosity, consistent with the idea that the NIR continuum emission mainly traces hot dust 
close to the dust sublimation radius, i.e. where the grains are heated above thei r sublimation 



temperature (T^uh ^ 1500 K) and d estroyed by the (stellar) radiation field (JTuthill et al. 



1/2 



20021 : iMonnier fc Millan-Gabetll2002l ). While the expected simple R^, oc Lj scaling law be- 



tween stellar luminosity L^ and NIR size i?* appears to hold throughout the low- to medium- 
luminosity part of the observed stellar sample, some very luminous early B-type stars exhib- 
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ited considerably smaller NIR sizes than predic te d by this relation (JMonnier fc Millan-Gabet 



2OO2I : Eisner et al.l 12004 : iMonnier et al 



2005h . iMonnier fc Millan-Gabetl J2002h suggested 



that this might be due to the presence of an inner gaseous disk, which shields the dust 
disk from the strong stellar ultraviolet (UV) radiation. Since this shielding would be most 
efficient for hot stars, it would allow the inner rim of the dust disk around B-type stars to 
exis t closer to the sta r . Several subsequen t studies favour "classica l" accretion disk mod- 



els (lEisner et al. 



2004 : IMonnier et al.l l2005l : IVinkovic fc Jurkid 120071 ) or "two-ring" models 



( JEisner et al.l 120071 ). in which the infrared emission contains contributions from the thermal 
emission of optically thick gas in the innermost disk regions. 



Following initial attempts by iHinz et al.l (|200ll ). who used single-dish nulling interfer- 
ometry to put upper limits of < 20 AU on the MIR size of some Herbig Ae stars, a large 
sample of HAeBes disks could be resolved with the VLTI/MIDI long-baseline interferometer 
at MIR wavelengths. In contrast to the NIR flux, which originates from hot dust in the 
innermost disk regions close to the dust sublimation radius, the MIR em ission also traces 
cooler (~ 200 — 300 K) dust further out in the disk. iLeinert et al.l (120041 ) determined char- 
acteristic dimensions of the 10 /im emitting regions for a sample of HAeBes, which ranged 
from 1 AU to 10 AU. 

Due to the quite limited uf-plane coverage of most existing infrared interferometric data 
sets, most published studies were only able to derive estimates of characteristic sizes and, in 
some cases, to look for possible elongation of the emitting region, but could not investigate 
the geometry of individual sources in detail. A more comprehensive ii iterferometric study o f 
one Herbig Ae star, HR 5999, was recently performed by our group (jPreibisch et al.ll2006l ). 
Based on a set of ten MIDI measurements at different projected baseline lengths and po- 
sition angles (PAs), modeling with 2-D frequency-dependent radiation transfer simulations 
provided relatively detailed information on the disk size and inclination. 



In this study, we combine, for the first time, NIR and MIR interferometric observations 
to constrain the spatial structure of the dust and gas environment around a Herbig Be star. 
Besides the wide spectral coverage (from ~ 2 to 13 yum), our interferometric data is also 
dispersed into several spectral channels (resolution R = A/AA ~ 30), allowing us to measure 
the wavelength-dependence of the visibility over the K- and A^-band and within spectral 
features such as the silicate emission feature around 10 /im. In combination with detailed 2-D 
radiative transfer modeling, this spectro-interferometric data set provides unique information 
about the inner circumstellar structures of this star. 



MWC 147 (alias HD 259431, BD+10 1172, HBC 529, V700Mon) is a well-studied Her- 
big Be star in Monoceros. There is some uncertainty concerning the distance and the 
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physical parameters o 



this star. From the analysis of the Hipparcos parallax data by 



van den Ancker et al.l (Il998l ). a lower limit on the distance of > 13 pc was derived, while a 
reanalysis suggested a distance of 290^84° pc (JBertout et al.lll999l ). This distance estimate, 



however, is in conflict with the apparent location of MWC 147 in the NGC 2247 dark cloud, 
which is a part of the cloud comp lex in the Monoceros OBI association at a distance of 
~ 800 - 900 pc Joiiver et alJll996l ). We assume a distance of 800 pc for MWC 147 (consis- 
tent with most other re cent studies) and use th e main stellar parameters as listed in Tab. [H 
which were taken from [Hernandez et al.l (120041 ). 



Numerous observational results strongly suggest the presence of a circumstellar disk 
around MWC 147. The object shows a strong infrared excess of about 6 mag at MIR wave- 



l ength s, clearly demonstrating the presence of circumstellar material. iHillenbrand et al. 
( 119921 ) fitted the SED of MWC 147 with a model assuming an accretion disk and estimated 
an accretion rate of Mace = 1-01 x 10^^ MqJt^^. MIR (10 /im and 18 /im) imaging ob- 
servations revealed an elongated diffuse emission component aroun d MWC 147 along PA 
50°, extending out to ~6" and contributing ~ 34% to the total flux (jPolomski et al.ll2002l ). 



Mannings! (1l994j ) determined the 1.1 mm flux of MWC 147 and estimated the mass in 
far disk/enve l ope t o be < 0.09 Mq. The study of the far-UV spectrum of 



the circumste 
MWC 147 by 



Bquret et al.l (120031 ) also suggested the presence of a flared circumstellar disk. 



Polomski et al.l (120021 ) imaged MWC 147 in the MIR and concluded that the star is sur- 
rounded by a mo derately flared disk and probably an extended envelope. Measurements by 
Jain et al.l (Il990l ) showed a significant amount of linear polarization (~1% along PA ~ 106°) 
but no wavelength-dep endence of the polariza tion. The high observed rotational velocity 
of vsini = 90 kms~^ (JBoehm fc Catalal Il995l ) suggests a high inclination of the rotation 
axis of MWC 147 with respect to the line-of-sight. This implies that the orientation of the 
circumstellar disk should be closer to edge-on than to face-on. 



Table 1. Stellar parameters for MWC 147 from iHernandez et al.l (120041 ). assuming 

Rv = 3.1. 



Parameter 



Value 



Spectral Type 




B6 


Effective Temperature 


Teff 


14125 K 


Luminosity 


L* 


155OL0 


Mass 


M* 


6.6 Mq 


Age 


t 


0.32 Myr 


Distance 


d 


800 pc 


Extinction 


Av 


1.2 mag 


Stellar Radius 


-R* 


6.63 Rq 
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Recently, iBrittain et al.l (120071 ) presented a high-resolution NIR spectrum of MWC 147, 
showing a strong Br7 emission hne. Using an empirica l relation between the Br7 luminosity 
and the accretion rate (as derived from UV veiling; Ivan den Anckerl 120051 ). they derive a 
mass accretion rate of Mace = 4.1 x 10""^ MQyr"^. As discussed in their Section 5.2, it is 
well possible that this method underestimates the true mass accretion rate. Adopting the 
stellar parameters used in our study (Tab. [1]) will also result in a larger value for the derived 
accretion rate. 

Evidence for a strong stellar wind from MWC 1 47 comes from the observed P Cygni 
profiles in several emission lines (JBouret et al.ll2003l ). A quantitative modeling of FUSE 
spectra revealed multiple absor ption components with different temperatures, consistent 
with a flared di sk interpretat i on (JBouret et al.ll2003l ). Based on the intensity ratio of infrared 
hydrogen lines, iNisini et al.l (119951 ) estimated a mass loss rate of 2.0 ± 0.4 x 10^^ MQyr~^, 
which is slightl y higher than the rn ass loss rates determined from radio observations (0.68 x 



10"^ Moyr-^ ISkinner et al.lll993h . 



The star has a faint visual compani on at a projected separ ation of 3'.'1 (~ 2 500 AU 



AR = 6.82, iBaineset al.ll2006h . While IVieira fc Cunhal (Il994[ ) classified MWC 147 as a 



spectroscopic binary with a period of about one year, this claim could not be confirmed in 
more recent observations (ICorporon &: Lagrangdll999l ). 



First interferometric measurements on MWC 1 47 were presen t ed bylMillan-Gabet et al. 



(j200ll ). providing an upper limit on the if -band size. lAkeson et al.l (120001 ) observed MWC 147 
with the Palomar Testbed Interferometer (PTI) and resolved its emission in the i^-band at 
baselines around 100 m. They derived a best-fit Gaussian FWHM diameter of 1.38 mas 
(=1.1 AU) in the K-hand. 



2. Observations and data reduction 



Details of all interferometric observations of MWC 147 used in this paper are summa- 
rized in Tab. [2l All visibility measurements were corrected for atmospheric and instrumental 
effects using calibrator stars observed during the same night. The calibrator stars as well 
as their assumed angular diameters are listed in Tab. [31 Fig. [1] shows the ■ut'-plane coverage 
obtained with these observations. 



Table 2. Observation log for interferometric observations on MWC 147. For more detailed information about the 

calibrator stars, we refer to Tab. [31 



Instrument 



Date HA Band/ Baseline Projected Baseline 

(UT) Spectral Mode Length [m] PA [°] 



Calibrators 



Ref. 



Near-Infrared 



IOTA/FLOUR 


1998 




H 


38 m 


~22 


~25 


see Ref. 




(1) 


IOTA/FLOUR 


1998 




K' 


38 m 


~21 


~ 15 


see Ref 




(1) 


PTI 


1999, 2000, 2003 


<0 


K 


NS 


105.1 


29 


HD 42807, 


HD 43042, 


(2), (3) 






>o 




NS 


98.9 


17 


HD 43587, 
HD 50692 


HD 46709, 




PTI 


2004 




K 


NW 


85.7 


76 


HD 43042, 
HD 50692 


HD 46709 


(3) 


PTI 


2003, 2004 


<0 


K 


sw 


78.6 


154 


HD 42807, 


HD 46709 








>0 




sw 


84.4 


143 


HD 50692 






VLTI/ AMBER 


2006-02-20 


03:14 


K/LR 


UT1-UT3 


101.0 


40 


HD 45415 






Mid-Infrared 


VLTI/MIDI 


2004-10-30 


08:49 


N/PRISM 


UT2-UT4 


89.4 


82 


HD 31421, 


HD 49161 




VLTI/MIDI 


2004-11-01 


05:23 


N/PRISM 


UT2-UT4 


55.9 


90 


HD 25604, 
HD 31421 


HD 49161, 




VLTI/MIDI 


2004-12-29 


06:42 


N/PRISM 


UT2-UT3 


46.5 


46 


HD 49161 






VLTI/MIDI 


2004-12-30 


02:33 


N/PRISM 


UT3-UT4 


59.6 


114 


HD 49161 






VLTI/MIDI 


2004-12-31 


04:26 


N/PRISM 


UT3-UT4 


61.5 


108 


HD 49161 






VLTI/MIDI (rej.) 


2005-01-01 


05:43 


N/PRISM 


UT3-UT4 


54.7 


106 


HD 31421, 


HD 49161 




VLTI/MIDI (rej.) 


2005-02-28 


00:04 


N/PRISM 


UT2-UT3 


39.3 


44 


HD 49161 






VLTI/MIDI 


2007-02-08 


05:41 


N/PRISM 


UT1-UT3 


102.0 


37 


HD 49161 






VLTI/MIDI 


2007-03-10 


03:12 


N/PRISM 


UT1-UT2 


56.4 


35 


HD 49161 







Oi 



References. — CP lMillan-Gabet et al]|200ll: (2) re-processing of data presented in lAkeson et al.ll2000l : (3) IWilkin fc Akesonir2003l . 
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Table 3. Calibrator star information for the interferometric observations presented in 

Tab.E 



Star 


K 


N 


Spectral 


I^UD.K 


liuD,iV 






[Jy] 


Type 


[mas] 


[mas] 


HD 42807 


4.85 


0.5 


G2V 


0.45 ± 0.03^ 




HD 43042 


4.13 


1.0 


F6V 


0.59±0.10'' 




HD 43587 


4.21 


1.0 


F9V 


0.48 ± 0.30<^ 




HD 45415 


3.02 


3.3 


G9ni 


1.06 ±0.02^ 




HD 46709 


2.62 


4.4 


K5ni 


1.66±0.02d 




HD 50692 


4.29 


0.8 


GOV 


0.56it0.10'' 




HD 25604 


2.03 


7.2 


KOIII 




1.96 ± 0.08° 


HD 31421 


1.41 


13.6 


K2IIIb 




2.58 ±0.15° 


HD 49161 


1.58 


7.2 


K4III 




2.88 ±0.17° 



Note. — The V-band magnitudes were taken from SIMBAD, the 
K-haad magnitudes from the 2M ASS point source catalog , and the 
Af-band (12 /xm) flux density from lHelou fc Walked lll988li . 

''UD diameter from lMalbet et al.l lll998l) . 

"U D diameter from the CHARM catalog IIRichichi &: PercheronI 
l2002ll . 

°UD diameter adopted from lPasinetti Fracassini et al.l (120011 ). us- 
ing the Hipparcos parallax of 51.76 mas measured for HD 43587. 

<*UD diameter from lMerand et al.l l l2005h . 

°UD diameter from the CHARM2 catalog JRichichi et al.ll2005h . 
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Fig. 1.— uv-plane coverage of the VLTI/MIDI, VLTI/AMBER, and archival PTI data. 
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Fig. 2. — Visibilities measured with MIDI as a function of wavelength. 
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2.1. VLTI/MIDI observations 



The MIDI interferometer (jPrzygodda et al.ll2003l : iLeinert et al.ll2004j ) at the ESO Very 
Large Telescope Interferometer (VLTI) records spectrally dispersed interferograms in the 
iV-band (8-13 /im). The MIDI observations of MWC 147 were carried out for ESO open 
time (OT) programmes 074.C-0181 and 078.C-0129 (P.I. Th. Preibisch), using the NaCl 
prism as dispersive element (providing a spectral resolution oi R = 30) and the HIGH-SENS 
instrument mode. In total, nine observations were carried out on five different baseline 
configurations, as listed in Tab. [2l 

We used the data reduction software package MIA+EWql| (Release 1.5.1) to extract 
visibilities from the MIDI data. This p ackage cont ains two independent data reduction 
programs, based o n a coherent (EWS, iJaffd |2004| ) and an incoherent approach (MIA, 
Leinert et al.ll2004j ). The reduction results obtained with both algorithms agree very well 
(within 3% for the calibrated visibility) for all data sets; with exception of the data sets 
from January/February 2005. Inspection of the acquisition image for these specific data sets 



iThe 



MIA-l-EWS 



software package 



available 



from 



the 



website 



http : //www . mpia-hd . mpg . de/MIDISOFT/ 



> 



0.8 



0.6 



0.4 



0.2 
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Fig. 3. — Wavelength-dependence of the visibihty measured with VLTI/AMBER and the 
visibihties measured with PTI using broad-band filters. 



revealed a poor beam overlap; thus, we rejected them from our further analysis. The inspec- 
tion of the acquisition images for the other data sets showed that the visual companion at 
a separation of 3'.'1 was not in the MIDI field-of-view (FOV) and, therefore, does not affect 
the measured visibilities. The wavelength-dependent calibrated visibilities of the remain- 
ing seven data sets are shown in Fig. [2J For the minir aum relative error o n the calibrated 
visibilities, we assume a conservative value of 10% (see iLeinert et al.ll2004j ). 



2.2. VLTI/AMBER observations 



AMBER flPetrov et al.ll2003l . 120071 ). the NIR beam-combiner of the VLTI, can combine 
the light from up to three telescopes, providing not only the visibility amplitudes but also the 
closure phase relation. The observations of MWC 147 were conducted within OT programme 
076.C-0138 (RI. Th. Preibisch) using the 8.2 m unit telescopes UT1-UT3-UT4. Spectrally 
dispersed interferograms were obtained in the low resolution (LR) mode {R = 35), which 
resolves the i^-band into 11 spectral channels. Due to problems with the fiber injection 
during that night, the flux reaching the AMBER beam combiner from UT4 was about a 
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factor of 3 lower than from the other telescopes. Therefore, clear fringes were detectable on 
only one of the three baselines (UT1-UT3), and no closure phase signal could be measured. 
Following the ESO service mode observation procedure, one data set of 5000 interferograms 
was recorded both on a calibrator star (HD 45415) and on MWC 147. The length and 
orientation of the projected baseline for this AMBER measurement (101 m, PA 40°) is 
similar to the measurement at the PTI-NS baseline, but adds information about the spectral 
dependence of the visibility within the ii"-band. 



The AMBER dat a were reduced wi th version 2.4 of the amdlido software employing 
the P2VM algorithm (JTatuUi et al.ll2006l ). Due to th e absence of a fringe tracker, a larg e 
fraction of the interferograms is of rather low contrast iPetrov et al.l (see discussion in 120071 ). 
Therefore, we removed those frames from our data set for which (a) the light injection from 
the contributing telescopes was unsatisfying; i.e., the intensity ratio between the photometric 
channels was larger than 4, (b) the atmospheric piston was larger than ^/^ of the coherence 
length A • i?, or (c) the fringe contrast was decreased due to instrumental jitter ( the 20% best 



interferograms were selected based on the Fringe SNR criteria, as defined in iTatuUi et al. 
20061 ). In Fig. [3] the calibrated K-hand visibilities derived from th e AMBER and PT I 



measurements are shown as a function of wavelength. As mentioned in lPetrov et al.l (120071 ). 
the accuracy of the absolute calibration of the visibilities measured with VLTI/AMBER is 
currently limited by vibrations induced by the UTs, which lead us to assume a minimum 
relative error of 3% for the calibrated visibilities. In contrast to the absolute ca libration, the 



wavelength-differential dependence of the visibility is insensitive to this effect ( IPetrov et al. 
2007h . 



2.3. PTI archive data 



MWC 147 wa s observed with the Palomar Testbed Interferom eter ('PTI. IColavita et al. 



19991 ) on the NS (lAkeson et al.l l2000l ) and NS & NW baselines (IWilkin fc AkesonI 120031 ) . 
Yet unpublished data for the SW baseline was retrieved from the PTI archive. To obtain a 
uniformly calibrated data set, we processed the new data set together with the previously 
published data using the V2Calib VI. 4 softwarqj. In the course of the cahbration procedure, 
we applied to the raw visibilities (which were estimated from the spatially filtered PTI 
spectrom eter output) the stan dard correction for coherence loss using the measured phase 
jitter (see IColavita et al.lll999l ). The individual PTI measurements were binned so that each 



^The amdlib software package is available from the website http://amber.obs.ujf-grenoble.fr 
•^The V2Calib software is available from the website http://msc.caltech.edu/software/V2calib/ 
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bin contains data sets covering less than 15° variation along the PA. Since the measurements 
on the PTI-NS baseline also cover a relatively wide range of PAs (~ 23°), we divided those 
measurements into two halves (depending on the PA) before averaging. As for the AMBER 
data, we assume 3% minimum relative error on the calibrated visibilities. 



2.4. Spitzer-IRS archive data 



In order to constrain the SED for our radiative transfer modeling as tightly as pos- 
sible, we obtained MIR spectra from the Spitzer Space Telescope Archive. These spectra 
were recorded on 2004-10-26 within GTO progr amme ID 3470 (P.I. J. Bouwman) using the 
Infrared Spectrograph (IRS, iHouck et al.l 120041 ). The data set consists of four exposures; 
two taken in the Short-High mode (SH, wavelength range from 9.9 to 19.6 /xm) and two in 
the Long-High mode (LH, 18.7 to 37.2 /im). Both modes provide a spectral resolution of 
R ~ 600. With slit sizes of 4'.'7 x 11'.'3 (SH mode) and ll'.'l x 22'.'3 (LH mode), IRS integrates 
flux from areas much larger than those collected in the spatially filtered MIDI spectrum. 
The spectra were pre-processed by the S13.2.0 pipeline version at the Spitzer Science Ce nter 
(SSC) and then extracted with the SMART software. Version 5.5.7 fJHigdon et al.ll2004l ). 



3. Results 



3.1. The MIR spectrum 
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Fig. 4. — Comparison of the measured MIDI spectrum and the spectrum extracted from 
archival Spitzer-IRS data. 
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In the overlapping wavelength regime between 10 and 13 fun, the MIDI and Spitzer-IKS 
spectra show good quantitative agreement, both in the absolute level of the continuum flux 
and in the spectral slope (see Fig. Hj). However, the IRS spectrum exhibits some line features 
which do not appear in the MIDI spectruno. As these emission lines are most pronounced at 
wavelengths of 11.0, 11.2, 12.8, 14.5, and 16.4 / xm, we attribute these features to the presenc e 
of Polycyclic Aromatic Hydrocarbons (PAHs, lAUamandola et al.lll985l : Ivan Dishoeckll2004l ). 
which were found towards a large v ariety of objects, including T-Tauri stars and HAeBe 
stars (jAcke fc van den Anckerll2004l ). For the strong and rather broad emission feature at 
11.2 fira, contributions from the 11.3 /im crystalline silicate feature are also possible. 

The prominence of the PAH lines in the Spitzer-IRS spectra with their 4'.'7 x 11'.'3 FOV 
and their absence in the MIDI spectrum with its much (~ 17x) smaller FOV suggests 
that the PAH emission comes predominantly from the outermost circumstellar environment 
and/or the surroun ding nebulosity of MWC 147, similar to w hat was found for other young 
stellar objects (e.g. van Boekel et al.ll2004bl . iRho et al.ll2006l ) or the outer regions of HAeBe 



disks ( jHabart et al. 



2004 



20061 ). 



The comparison of the Spitzer-IRS spectrum with the fluxes measured by IRAS (JHelou fc Walker 



19881 ) showed that the IRAS fluxes are systematically higher. This is likely related to the 
larger beam size of IRAS, which includes signiflcant amounts of emission from the ambient 
NGC2247 nebula. For our modeling in Sect. HI we will therefore treat the IRAS fluxes as 
upper limits only. 



3.2. The correlated MIR spectrum — indications of grain growth 



The visibilities measured with MIDI show significant variations along the recorded wave- 
length range. In particular, we detect a drop of visibility within the 10 /xm silicate feature. 
A similar behav i or ha s alre ady been observed for se veral HAeBe stars; e.g., in the samples 
of iLeinert et all (|2004J ) and Ivan Boekel et all (l2004al ). 



As the silicate emission featu re is generally attributed to the presence of rather small 
silicate grains (r < 0.1 /im, see e.g. Ivan Boekel et al.ll2004al ). it is possible to probe the radial 
dust mineralogy by comparing the correlated spectrum at various baseline lengths with the 
total spectrum Ftot- The correlated spectrum Fcorr corresponds to the flux integrated over the 



"^In order to validate that the non-detection of the strong emission hne around 11.2 /im in the MIDI 
spectrum is not only due to MIDI's lower spectral resolution, we convolved the Spitzer-IRS spectrum to 
the resolution of the MIDI PRISM. Since the line is still clearly visible in the convolved spectrum, we are 
confident that no significant line emission is present in the spatially filtered MIDI spectrum. 
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Fig. 5. — Total MIDI spectrum and the correlated flux at the 55.9 m (corresponding to a 
spatial scale of < 35 AU) and 89.4 m (< 20 AU) baselines. 



spatial area unresolved by the interferometer for a particular baseline length B. Therefore, 
for each baseline length B, the correlated flux Fcott{B) can be computed by multiplying 
the total spectrum measured by MIDI in the photometry files with the visibility measured 
for a certain baseline. In order to probe the radial dependence of the dust mineralogy, 
measurements taken at similar PAs should be used in order to avoid contaminations by 
changes in the source geometry. We therefore choose the measurements from 2004-11-01 
{B = 55.9 m, PA=82°) and 2004-10-30 {B = 89.4 m, PA=90°), as they have very similar 
PAs. The comparison of the correlated spectra for these baselines with the total spectrum 
(Fig. [5]) shows that the 10 /xm silicate feature flattens out with increasing resolution. This 
change in the correlated spectrum might indicate spatial variations in the dust co mposition, 



with the more evolved dust grains (i.e. larger grains with a weaker silicate feature; iMin et al. 
20061 ) in the innermost disk regions. 



3.3. Geometric model fits 

Since the imaging capabilities of the current generation of infrared interferometers are 
rather limited, the measured interferometric observables are often used to constrain the 



-16- 



parameters of a model for the object morphology. In most studies presented until now, 
either purely geometric profiles (in particular uniform disk (UD) and Gaussian profiles) or 
physically motivated geometries such as ring profiles or analytic accretion disk models with 
a power law temperature distribution were employed. Ring models are justified by the 
theoretical expectation that most of the NIR emission origiriates from a rather small region 



aroun d the dust sublimation radius (e.g., iMillan-Gabet et al.ll200ll : lMonnier fc Millan-Gabet 
2002h . 



A common problem in applying simple geometric models is that the observed emission 
does not originate exclusively from the circumstellar material: a certain fraction comes 
directly from the central star and contributes a spatially (nearly) unresolved component, 
and the existence of extended background emission, which is fully resolved, is also possible. 
For the model fits, one therefore has to specify which fraction of the total flux Ftot at 
any wavelength has to be attributed to the different spatial components. The stellar flux 
contribution /star/tot (A) = -Fstar/-Ptot is often estimated from the SED, while the extended 
component /ext/tot(A) = -FcxtZ-^tot is usually assumed to be zero. These assumptions are, 
however, associated with a considerable uncertainty. 

To allow c omparison with earlier NIR interferometric studies on MWC 147, we keep the 
flux ratios from lMillan-Gabet et al.l(l200ll ). namely /star /tnt (2. 1 p,m) = 0.16, and fp:yt/tct ('2.l /im) 



0.0 for the aii a lytic fits. The same flux ratio was assumed by lWilkin fc AkesonI (120031 ) . while 
Akeson et al.l (120001 ) used /star/tot(2.1 /xm) = 0.10. At MIR wavelengths, the stellar contri- 
bution is likely to be negligible; i.e., /star/tot(10 /xm) ?» 0. This can be concluded from the 
SED shown in Fig. [6l where the infrared excess exceeds the stellar flux by a factor of ~ 280 
at 10 fj,m. 



3.3.1. Characteristic source size and elongation 

To obtain a flrst estimate of the object size, we flt the most common analytic proflles 
to our interferometric data: Gaussian, UD, and ring proflles. For mathematical descr i ptions 



of these proflles, we refer to iKraus et al.l (120051 . UD proflle) and iMillan-Gabet et al.l (12001 



Gaussian & ring proflle). We assume un iform bright rings wit h an average diameter and a 
flxed width of 20% to be consistent with lMonnier et al.l (120051 ). As the apparent object size is 
expected to change with wavelength, we fltted these proflles to subsets of our data, covering 
wavelength ranges around 2.2 yum, 8.6 /im, 11.0 yum, and 12.5 /xm. The visibilities measured 
in these sub-bands were fltted to the visibility proflles using a Levenberg-Marquardt least- 
square fltting algorithm, taking the chromatic change in resolution within the bandwidth 
into account. 
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Fig. 6. — SED for MWC 147 containing the measured averaged, weighted MIDI spectrum, 
the archival Spitzer-IKS spectrum and data from literature. Comparing the de-reddened 
SED with the Kurucz stellar atmosphere model for a B6-star reveals excess emission both 
in the UV and IR wavelength regime. The UV excess emission migh t originate in mag- 
netospheric accretion (see the general models by iMuzeroUe et al.l 120041 ). while the infrared 
excess likely indicates the presence of circumstellar du st and optically t h ick ga s (see mod- 
eling in Sect. H]). We included photometric dat a frora [Thompson et al.l (Il978l . 156.5 nni, 
196.5 nm 236.5 nm, 274.0 nm), |Wesselius_elalJ JlQsd 150 nm, 180 nm, 220 nm, 250 nm 
330 nm). [Jnienbrand et aP Jl992l. U . B, V, RC, IC^JJLJC L, M N, Q-bands) . Furon et al 



J1993I. l^-ba nd \. fegret et al.l Jl992l. 5-band), Ihos et a 



Skrutskie et all (l2006l . 2MASS, J, H, /sT.-band) , [Cohenl (119731 . 2.2 /im, 3. 5 /im, 3.65 um 



tod, Vt, Bt, Rj, /j-band). 



.8 jUm , 4.9 /im, 8.4 /im, 8.6 /im, 10.8 /im, 11.0 /im, 11.3 /zm, 12.8 /zm, 18 /im) . iPolomski et al. 



2002 



4.74 Mm 7.91 /im, 8.81 /im, 10.27 /im, 11.70 /im, 12.49 /im, iV-ban d, 18.17 ^m), 



Egan et all fll999l . MSX 2.6, 4.29 /im, 8.28 /im, 12.13 /im, 14.65 /xm, 21.34 ^m), 



fll987l . 2.2 /im, 3.85 /zm, 8.65 / xm, 9.97 ^m 10.99 /zm, 11.55 a mi . iHelou fc Walker 



12 /im, 25 /im, 60 /im, 100 /im), [Caseyl fll99ll . 160 /xm, 370 /im) . iManningsl fll994l . 450 /im) 



Berril 



i et al. 



(11988 
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Fig. 7. — Polar diagram showing the best-fit circular and elliptical geometries. To derive 
the characteristic object size along different PAs (dots with errorbars), we fitted ring profiles 
to the PTI i^-band visibilities (left panel) and to MIDI visibilities (right panel) which were 
averaged over five spectral channels around the silicate feature (10.3 — 11.7 /xm) and in the 
surrounding continuum ranges 8.1 — 9.1 /xm and 11.9 — 13.1 /xm. 



-19- 



18 
16 
14 
12 



en 
re 

£ 10 

i 8 
re 

6 
4 



2 - 




2.1 2.2 



2.3 2.4 



.,../- 



T °= r~^'^ disk (Normalized to NIR) 
T °c r"'^^'^ disk (Normalized to NIR) 
T oc r"^'^ disk (Normalized to MIR) 
T oc r disk (Normalized to MIR) 



t? 



„ 




6 8 

Wavelength [nm] 



10 



12 



14 



Fig. 8. — Wavelength-dependence of the measured characteristic size over the H-, K- and 
A'^-band, including IOTA, PTI, AMBER and MIDI measurements. From our seven MIDI 
measurements, we choose the data set from 2007-02-08, since this measurement was taken at 
very similar baseline length and PA (102.0 m/35°) as the AMBER data set (101.0 m/40°), 
which makes these data sets particularly well suited to study the radial disk structure without 
contamination by the detailed source geometry. The IOTA if -band measurement provides 
only an upper size limit. To compute the characteristic size, Gaussian intensity profiles were 
assumed. For comparison, we show the w avelength-depende n t size corresponding to the 
commonly applied analytic disk model from iHillenbrand et al.l (Il992l ). We scaled these disk 
models to match the measured NIR (dashed curves) or MIR size (dotted curves). In all cases, 
it is evident that these analytic models cannot describe the measured wavelength-dependent 
size well. 
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The fits were performed for tlie case of circular symmetry (e.g. a disk seen face-on) and 
for elliptical structures (e.g. an inclined disk). The obtained diameters and goodness-of-fit 
values {xl = S [(^^ ~ ^m)/^v] l^v, with V"^ being the measured squared visibility, V!^ 
the squared visibility computed from the model, and Ny the number of measurements) are 
given in Tab. H] and [51 

The Xr values already indicate that the elliptical geometries are a better representation of 
our data than the circular models. In order to illustrate this object elongation, we show the 
corresponding geometries in Fig. [71 Since the detection of object elongation requires strict 
uniformity in the observational methodology, we did not include the single-baseline AMBER 
measurement because the mixture of broadband and spectrally dispersed interferometric 
observations might easily introduce artefacts. 

While the elongation is only marginally evident in the PTI NIR measurements, it is more 
significant in the MIR data, although the limited PA-coverage of the available MIDI data 
(covering ~ 80° in PA) prevents us from measuring the precise axis ratio. The effect that the 
deviations from circular geometry seem to be stronger at MIR wavelengths might indicate 
that for the NIR model fits, the assumed stellar flux (/star /tot (2.1 /xm) = 0.16) underestimates 
the real flux contributions from a spatially unresolved region (in agreement with our radiative 
transfer modeling results in Sect. 14.41) . For the MIR-interferometric data, it is interesting 
to note that the elongation fou nd agrees well w i th th e one seen in the 11.7/18.2 /xm color 



temperature map published by [Polomski et al.l (120021 ). Although these color temperature 



maps show structures on scales of several arcseconds (e.g., on scales a hundred times larger 
than our interferometric data), their orientation (~ 50°) and rough axis ratio are similar to 
those of the structure seen in our MIDI observations. 

Table 4. Model fits for circular geometries. 

_R'-band 9 /im 11 ^in 12.5 /xm 



Diameter Xr Diameter 

[mas] [mas] 



Uniform Disk (UD) 2.6 ± 0.2 1.44 12.8 ± 2.0 1.38 17.3 ± 2.2 0.95 18.1 ±2.2 0.89 

Gaussian 1.6 ±0.1 1.30 7.9 ± 1.4 1.19 10.8 ± 1.5 0.75 11.2 ± 1.5 0.90 

Ring 1.8 ±0.1 1.51 8.7 ±1.4 1.49 11.8 ±0.3 1.07 12.4 ±1.3 0.90 



Note. — For the ii'-band fits, we attribute 16% of the total flux to the unresolved stellar component. At a 
distance of 800 pc, 1 mas corresponds to 0.8 AU. 



Table 5. Model fits for inclined geometries. 
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Note. — The angle i denotes the inclination computed from the fitted axis ratio, assuming an underlying circular source structure. However, due to the 
unknown flux contribution from the unresolved central region (i.e. stellar flux, for which we assume /star/tot (NIR) = 0.16 and /star/tot (^''IR) = 0-0, and maybe 
additional accretion luminosity), and the incomplete PA-coverage of the MIDI data set, we consider this value rather uncertain. 
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3.3.2. Wavelength- dependent characteristic size 

Our combined MIDI/AMBER/PTI/IOTA data set also allows us to study the wavelength- 
dependence of the apparent size. For this, we fitted the visibility measurement in each indi- 
vidual spectral channel with the analytic formula for Gaussian intensity profiles (the result 
does not depend strongly on the assumed profile). The determined diameters are shown in 
Fig. [HI The increase of the apparent size with wavelength is usually interpreted as a conse- 
quence of the radial temperature profile for the circumstellar material (i.e. material at larger 
distances from the star is cooler). However, as will be quahtatively discussed in Sect. 13. 4[ 
the increase of the characteristic size with wavelength is much steeper than expected for sin- 
gle power-law temperature accretion disks with an inner hole, suggesting that these models 
cannot successfully explain the observed wavelength-dependence of the emitting structure. 
In the context of our radiative transfer modeling in Sect. HI we will discuss that this dis- 
crepancy can be explained assuming the presence of an additional, more compact emitting 
region which dominates the emission at NIR wavelengths, but whose contribution becomes 
less important at MIR wavelengths. 



Akeson et al.l (120001 ) examined the possibility whether the measured visibilities could 



indicate the presence of a close companion. Since PTI observed MWC 147 several times over 
a period of ~ 4 yrs and found no significant variations on the NS-baseline, we consider the 
binary scenario as very unlikely. For example, assuming, just for the sake of argument, a 
total system mass of 7 Mq and a semi-major axis of 10 AU, this would give an orbital period 
of ~ 12 yrs and should result in significant visibility variations over the covered 4 yrs. 



3.4. Comparison with analytic disk models 



Analytical models, both of passiv e irradiating circumstellar di sks (JFriedjungl Il985l ) as 
well as viscous, actively accreting disks (JLynden-Bell fc Pringlelll974j ). predict that the radial 
temperature profile of YSO disks sh ould follow a simple power-law T(r) oc r~". Most studies 
infer a power law in dex of a = 3/4 (JMillan-Gabet et al.ll200ll : lEisner et al.ll2005l ) or 1/2 (e.g. 
Leinert et al.ll2004l ). Using this radial temperature power law and the assumption that each 
disk annulus radiates as a blackbody, we can compute the wavelength-dependence of the disk 
size corresponding to a certain analytic model. To model the wavelength-dependence of the 
size of an accretion disk model with constant power law, we simulated disks with T{r) = r~", 
where r is chosen to be in units of the disk inner radius with a temperature at the inner dust 
rim of T(l) = 1500 K. The outer disk radius is chosen so that the temperature drops well 
below 100 K. After computing the intensity profile, we determ ine the disk si z e ajX ) for these 
analytic disk models using the half-light radius definition by iLeinert et al.l (120041 ). Finally, 
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a) AnalTtic Disk Model (HUleiibiaml et al. 1992> 

Diigt + Gas 

b) VEltT-ltlM: Flared disk with vertical rim 



^H Dust + Gas 




Gas 
(ojil. iTiiek »r iTiiri) 



c) CLRV-RDl: Flared disk tvitli curved ptiffed-up rim 





Fig. 9. — Illustration of the co nsidered disk mod e l geom etries: (a) "Classical" accretion disk 
geometry as considered in the [Hillenbrand et al.l (Il992l ) model (see Sect. 13. 4p . incorporating 
an optically thick gas & dust disk, (b) In our models VERT-RIM and VERT-RIM-ACC, the 
dust component of the outer flared dust disk is truncated at the dust sublimation radius, 
while the gas located at smaller radii is either optically thin (Model VERT-RIM, Sect. 14.31) 
or optically thick (Model VERT-RIM-ACC, Sect. KM- (c) Due to direct stellar heating, the 
scale height at the inner dust rim might be increased, resulting in a puffed-up inner rim. 
Furthermore, various effects, such as a pressure- dependent dust evaporation temperature 
and dust segregation towards the disk midplane, will cause the rim to curve. We investigate 
the influence of such a curved rim shape in our models CURV-RIM (Sect. 14. 3p and CURV- 
RIM-ACC (Sect.lOj). 
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we scale the half-light diameter to fit the NIR (or MIR) size measured on MWC 147 and 
show the resulting a(A)-curves for the afore-mentioned representative values for a (Fig. [8]). 



Hillenbrand et al.l (Il992l ) fitted the SEDs of HAeBe stars with analytic disk models 
that combined the temperature profile of the reprocessing and actively accreting disk com- 
ponent. This class of "classical" optically thick, geometrically thin accretion disk models 
(Fig. [9]a) was also favoure d by various authors to model the SEP and NIR broadband vis i- 
bilities of Herbig Be stars (JEisner et al.ll2004l : iMonnier et al.ll2005l : IVinkovic fc Jurkig 120071 ). 
Vinkovic fc Jurkid (120071 ) modeled a large set of archival interferometric data and concluded 
that such classical accretion disk models are consistent with the observed visibilities of Her- 
big Be stars. Therefore, we examined whether such models can reproduce our combined 
NIR/MIR spectro-interferometric data on MWC 147. In order to fit this analytic model 
with precisely the same procedure as for our radiative transf er modeling (Sect . 14. 1.21) . we 
computed the temperature profile using the equations given in iHillenbrand et al.l (119921 ) and 
included the blackbody emission from an infinitely thin disk in our radiative transfer grid. 
Using the stellar paramet ers given in Tab. [ T , we found that a disk model with the same 
parameters as derived by [Hillenbrand et al.l (1l992l ) for this star (i.e. inner and outer disk 
radii of -Rj„ = 12 R^ and Rout = 62 AU) provides a good fit to the observed SED for an 
intermediate disk inclination angle (45°, see Sect. I3.3.ip and assuming a mass accretion rate 
of Mace ~ 4 X 10~^ M0yr~^. As shown in Fig. [TOl this model can also reproduce the ab- 
solute level of the NIR visibilities, but it fails to reproduce the spectral dependence of the 
VLTI/AMBER visibilities and also underestimates the MIR size, resulting in a rather high 
X^ of 5.56. 

It is obvious that these analytic models cannot reproduce the measured NIR and MIR- 
sizes simultaneously. To understand the strong increase in the apparent size with wavelength, 
contributions from the following effects might be of importance: 

(a) Analytic disk models generally do not include the effects of scattered light, which 
can provide significant heating of the outer parts of the disk and, thus, increase the apparent 
disk size at MIR wavelengths. 

(h) It has been suggested that the disks around YSOs may not be fiat but fiare with 
increa sing radius. Such a fiaring i s expected from vertical hydrostatic equilibrium considera- 
tions (JKenyon fc Hartmannlll987l ). As a result, the outer disk regions intercept more stellar 
fiux, resulting in an increased luminosity and apparent size at MIR wavelengths. 

(c) The fiux contribution from an extended cold envelope /onv/tot('^) might be non- 
negligible in the MIR. 



(d) The NIR size might be underestimated, if the amount of NIR emission originating 
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from close to the star is inadequately estimated (e.g. due to a biased /star/disk (A) or additional 
accretion luminosity). 

This enumeration illustrates that the currently routinely applied analytic disk models 
contain several problematic points. A more physical and consistent approach requires de- 
tailed radiative transfer modeling. In the next section, we will present 2-D radiative transfer 
modehng for MWC 147. 



4. 2-D Radiative Transfer Simulations 
4.1. Modeling procedure 

4-1.1. Monte Carlo radiative transfer simulations 

For a physical modeling of our interferometric data, we employ the radiative transfer 
code mcsimjmpi (author: K. Ohnaka), which solves the radiative t ransfer proble m self- 
consistently using a Monte Carlo approach. This code is described in 
and has been ap p hed f or the interpretation of interferometric data in 



Ohnaka et al. 



Ohnaka et al. 



(2006) 

jiood) 



and iHonig et al.l (l2006l ). In short, the stellar flux is treated as a finite number of photon 
packets which are emitted in arbitrary directions. While propagating through the cells of 
the simulation grid, the photon packet can be either scattered or absorbed. The probability 
of these events is given by the density and the optical properties of the dust in each par- 
ticular cell. For scattering events, the propagation direction of the photon packet changes 
anisotropically, according to the Henyey-Greenstein function (this allows us to also properly 
treat scattering on large grains). For each absorption event, energy is deposited into the 
cell while the packet is isotropica lly re-emitted immediately . The temperature of the cell 



is corrected using the scheme by iBjorkman fc WoodI (J200ll ). resulting in a self-consistent 



determination of the dust temperature distribution. After tracing the propagation of a large 
number of photon packets through the simulation grid, the SED is computed by summing 
the flux from all packets. The code is parallelized using the LAM/MPI library, which allows 
the user to distribute the Monte Carlo computation on a large number of computers within 
a network. We extended the original mcsimjmpi code by adding an option to include the 
thermal emission from optically thick gas (see Sect. 14.41) . 



As we require particularly high spatial resolution in the inner disk region to properly 
resolve the structure of the inner dust rim (at scales of a few AU) but also need to include 
structures with large radial extension (10 000 AU scale), we employed a spherical grid with 
logarithmic radial grid spacing. The number of radial cells was chosen to be 500, while 
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Fig. 10. — Analytic accretion disk model, assuming the [Hillenbrand et al.l (119921 ) disk tem- 
perature profile {Xr = 5.56). a) and b) show the corresponding brightness distribution for 
two representative NIR (2.25 ^m) and MIR (10.0 ^m) wavelengths, c) shows the SED for 
various inclination angles, whereas d) gives the SED for the best-fit inclination angle and 
separates the flux which originates in stellar photospheric emission and from the accretion 
disk. Finally, e) and /) depict the NIR and MIR visibilities. 
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Table 6. Parameters and fitting results for our 2-D radiative transfer models. 



SHELL 

(Fig. HI} 



VERT-RIM 

(Fig. m 



CURV-RIM 
(Fig.O 



VERT-RIM-ACC 
(Fig.EJ 



CURV-RIM-ACC 
{Fig.linil 



Outer radius [AU] 

Dust density at 10 AU^ [g-cm~' 

Grain size distribution 

Radial power law exponent p 

Vertical power law exponent q 

Mass accretion rate"^ [Mqjt" 



100 


100 




100 


100 


100 


2.8 X 10-21 


1.1 X 10- 


-18 


1.6 X 10-18 


9.4 X 10-19 


8.0 X 10-19 


Small 


Small & Large 


Small & Large 


Small & Large 


Small & Large 


3/2 


15/8 




15/8 


15/8 


15/8 


- 


9/8 




9/8 


9/8 


9/8 


- 


- 




- 


7.0 X 10"" 


6.5 X 10-6 


- 


50° 




40° 


50° 


40° 


- 


50° 




30° 


80° 


110° 


316.7 


135.2 




119.4 


2.07 


2.71 


33.8 


3.72 




1.60 


0.79 


0.95 



Inclination 
Best fit PA 
xl NIR 
X?MIR 



xi total 



80.2 



25.3 



20.9 



0.99 



1.24 



Note. — All models include an extended outer envelope in order to fit the MIR to FIR BED (see Sect. I4.1.2| l. The envelope is 
composed of small dust grains, with a density distribution p oc r-i" with po = 3 X 10— ^^ g-cm- ^ at ro = 10 AU (see equation [2]| . The 
outer cutoff radius r^utoff = 40 000 AU was chosen such that the dust temperature drops below 10 K. 
For all models, we assume a minimum dust density of 10-^^ g-cm-'^. 

''Dust Density, as measured at a radius of 10 AU in the disk midplane. 

I'Dust Composition — Mixture of warm silicates l lOssenkopf et al.lll992r) and amorphi 
size distribution follows n{a) oc a 
ciniin = 1.0 pm and amax = 1 000 pm 



^■^. For Small Grains we choose On 



ous carbon lIHa 
0.005 pm and amax = 1.0 pm, whereas for Large Grains we use 



epliii) (to equal parts). The 



'^Inner Gaseous Disk — The inner gaseous disk component is modeled to be geometrically thin, optically thick, and to extend from the 
co-rotation radius iJcorot to the dust sublimation radius (see Sect. l4?4t . 
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Fig. 11. — Radiative transfer model computed for MWC 147 assuming a Spherical Shell 
geometry (Model SHELL, Xr = 80.2, see Sect. 14.21 for a model description). For the 
model parameters, see Tab. El In a) and b), we show the dust density and the temperature 
distribution, c) and d) show the ray-traced images for two representative NIR (2.25 /xm) 
and MIR (10.0 /xm) wavelengths, e) shows the SED for various inclination angles, whereas 
/) gives the SED for the best-fit inclination angle and separates the flux which originates in 
stellar photospheric emission, thermal emission, dust irradiation, and accretion luminosity. 
Finally, g) and h) depict the NIR and MIR visibilities computed from our radiative transfer 
models. 
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Fig. 12. — Similar to Fig. [TT], but showing the radiative transfer model for a Flared Ke- 
plerian disk geometry with a vertical inner rim (Model VERT-RIM, Xr — 25.3, see 
Sect. 14.31 for a model description). 
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Fig. 13. — Similar to Fig. [TT], but showing the radiative transfer model for a Flared Keple- 
rian disk geometry with a curved puffed-up inner rim (Model CURV-RIM, Xr = 20.9, 
see Sect. 14.31 for a model description). 
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Fig. 14. — Similar to Fig. [TTl but showing the radiative transfer model for a Flared Kep- 
lerian disk geometry with a vertical inner rim, including accretion (Model VERT- 
RIM- ACC, Xr = 0.99, see Sect. 141 for a model description). 
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Fig. 15. — Similar to Fig. [TT], but showing the radiative transfer model for a Flared Keple- 
rian disk geometry with a curved pufFed-up inner rim, including accretion (Model 
CURV-RIM-ACC, Xr = 1-24, see Sect. lU for a model description). 
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the latitudinal grid resolution is 1°. We used 5 x 10^ photon packets per simulation, which 
ensures sufficient statistics for Monte Carlo radiative transfer. 



As the input stellar spectrum, we use a Kurucz stellar atmosphere model (JKurucallOTOl ) 
for a B6-type star of solar metallicity (Teg = 14 00 K, logq = +4.04). For the optical 



dust pr operties we us e a mixture of warm silicate (lOssenkopf et al.l Il992l ) and amorphous 
car bon (Hannerl 119881) grains. The grain size distribution follows the dependence suggested 
by iMathis et al.l (Il977l . n{a) oc a~^-^, where a^^^ < a < CLmax)- For the outer envelope, 
we use small grains (amin = 0.005 /im; Omax = 1-0 /^m), whereas for the inner disk region, 
we mix in larger dust grains (amin = 1-0 yum; Omax = 1000 /xm). The species of large and 
small grains are treated separately in the radiative transfer computations. A minimum dust 
density of 10"^^ g-cm~^ is assumed. 

In order to avoid unrealistically sharp cutoffs at the outer edges of the considered density 
distributions, we used a Fermi-type function -F(rcutoff) to obtain a smooth truncation of the 
density distribution around the radius rcutofr: 

-1 



FiVcntos) 



1 + exp 



r - rcutoff \ 



e • '"cutoff / . 
where e = 0.05 defines the relative width of the transition region. 



(1) 



4-1.2. Joint fitting of the SED and wavelength- dependent visibilities 

Once the radiative transfer computation is completed, a ray-trace program is used to 
compute the SED and synthetic images for any wavelength of interest. For this project, we 
compute synthetic images at 3 wavelengths in the ilf-band, 3 wavelengths in the if-band, 
and 8 wavelengths in the A^-band. Finally, visibilities are computed from the simulated 
images for the points of the ww-plane covered by the data. In order to treat the visibility 
slope in the spectro-interferometric observations properly, we compute the visibility for each 
spectral channel of the MIDI and AMBER observations separately, using the synthetic image 
computed for a wavelength as close as possible to the central wavelength of the spectral bin. 

Besides the interferometric observables, we use the SED of MWC 147 (as shown in 
FigThe. [6l including the Spitzer-IR^ spectrum and photometric data from the lit erature) to 
constrain our radiative transfer models. The contribution of the visual companion (JBaines et al, 



20061 . separation 3'.'1) to the SED is unknown. Furthermore, the photorn et ric data may be 



contam inated by the ambient reflection nebula NGC 2247 (jCaseyl Il99ll ). iPolomski et al. 



(I2OO2I ) studied the appearance of the circumstellar environment of MWC 147 at MIR wave- 
lengths and found an extended structure of ~ 12" diameter (~ 9 600 AU) at 10 yum. 
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In the course of our modeling of the SED of MWC 147, we found that disk models are 
not able to reproduce the observed SED long- ward of 15 yum (see Fig. [6]). The large observed 
fluxes at mid- and far-infrared (FIR) wavelengths require the presence of an extended enve- 
lope. We tried various geometries for this envelo pe, including ro t ating , infalling envelopes 



with and without polar outflow cavities, as used by lWhitney et al.l (120031 ). However, the best 



agreement was obtained for a simple power-law dust distribution with a rather fiat slope 

Penv(r) = PO ( — 1 X F(rcutoflf), (2) 

where po is the dust density at some characteristic radius tq. 

The fitting of the SED and the visibilities was done in the following way: First, we fixed 
the geometry of the extended envelope to match the MIR to FIR SED. Then, we varied the 
geometry for the inner component (Sect. 14.21 to I^l3l) to obtain the best fit to the NIR and 
MIR SED. The SED could be reproduced well with any of the geometrical models considered 
below, again demonstrating the highly ambiguous nature of SED fits. Significant deviations 
between the model SEDs and the observed SED are only apparent in the UV spectral range, 
between ~ 150 nm and ~ 300 nm. The observed UV excess is related to accretion activity 
and originates probably from the shocks where accreted material crashes onto the stellar 
surface. As these accretion shocks are not included in our modeling, we ignore the SED 
deviations in the UV spectral range. 

Finally, the agreement of the model with the interferometric observables was tested. 
The analytic model fits presented in Tab. [S] indicate that the disk should be notably inclined 
(~ 40 to 60°); therefore, we computed each radiative transfer model for different inclinations 
(from edge-on to face-on with an interval of 15°) in order to find the best agreement both 
with the spectral shape of the SED and the visibilities. In total, about one thousand radiative 
transfer models have been computed to identify the best-fit models presented in Sect. l4.2l to 



4.2. Model SHELL: Spherical shell geometry 



Miroshnichenko et al.l (119971 ) proposed that optically thin shells can reproduce the SED 



of HAeBe stars. For our spherical shell model we assumed a density distribution given by 



x ^P 



Pshell(^) = PO — X ^(^cutoff)- (3) 

with p = 3/2. As can be seen in Fig. [TT|, the shell model can reproduce the measured SED 
quite well, but completely fails to reproduce the NIR and MIR visibilities, both the absolute 
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level of the visibilities as well as their spectral dependence. Due to this very poor fit to the 
interferometric observables [x^ = 80.2), we can reject this geometry. In addition, spherical 
models cannot reproduce the elongation revealed by the analytic model fits (see Sect. 13.31) . 



4.3. Models VERT-RIM & CURV-RIM: Passive dust disks 



One solution for a passive disk which is gravitationally dominated by the central star, 
vertically isothermal, and in vertical hydrostatic equilibrium is given by a Keplerian disk 
geometry (jShakura fc Syunyaevlll973l ) with a density distribution 



Pdisk{r,z) 
hz{r) 



Po 



To 



-p 



exp 



TT / Z 



X F(rcutoff) 



horo 



ro 



(4) 
(5) 



where r is the radial distance in the mid-plane, z is the height above the disk mid-plane, and 
hz is the vertical pressure scale height. Hq is the relative geometrical thickness of the disk 
at the characteristic radius Tq. To reproduce t he SEP of MWC 147, we fi nd best agreement 
with a slightly flared disk geometry [q = 9/8, iKenyon fc Hartmannlll987l ) and with a radial 
density power law exponent p = 15/8 (corresponding to a surface density law S(r) oc r"'^/'^). 

In the model VERT-RIM, the inner disk edge of the density distribution is simply trun 



cated, resu lting in a vertica l wall at the dust sublimation radius (Fig. ^). iNatta et al. 



(I2OOII ) and iDuUemond et al.l (I2OOII ) proposed a modification of the disk geometry; namely, 
a "puffed-up" inner rim. As the dust at the inner rim (at the dust sublimation radius) is 
directly ex posed to the st e llar ra diation, the scale height in this region will be significantly 
increased. Ilsella fc Nattal (120051 ) pointed out that Tgubi depends on the gas density, which 



results in a higher sublimation temperature in the disk midplane than in the disk atmo- 
sphere, causing a curved rim shape. Recently, iTannirkulam et al.l (120071 ) investigated also 
the effect of dust segregation on the rim shape, considering two dust species with different 
grain sizes which are distributed over different scale heights. Both the pressure-dependent 
dust sublimation temperature, as well as the dust settling effects, result in a curved shape of 
the inner dust rim. In order to investigate the in fluence of this curved rim on our model flt, 
we use the analytic dust segregation model from ITannirkulam et al.l (120071 ) to compute the 
dust scale height hz at the evaporation front as a function of distance along the disk mid- 
plane. We use the Planck mean opacities corresponding to the two dust species descripted 
in Sect. 14.1.11 and assume th at larger grain species has settled to 60% of the scale height of 
the smaller grain population (ITannirkulam et al.ll2007l ). The computed scale height for the 
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curved rim is combined with the scale height of the outer flared disk, providing the input 
density distribution for our radiative transfer simulation of a disk model with curved inner 
rim (Model CURV-RIM, Fig.Elc). 

For all models with irradiated disks, we find that especially the NIR model visibilities are 
always much smaller than the measured visibilities (see Fig. [T2]and [T3l) . A similar, although 
less pronounced, deviation was found in the MIR visibilities. Therefore, the radiative transfer 
modeling confirms and quantifies the general tendency already observed in the geometric 
model fits; namely, that without considering accretion luminosity, the measured NIR radius 
of ~0.7 AU (see Tab. H]) is a factor of 4 smaller than the dust sublimation radius. We conclude 
that although passive irradiated circumstellar disk models are able to reproduce the SED 
of MWC 147, these models are in strong conflict with the interferometric measurements 
(resulting in Xr = 25.3 for Model VERT-RIM and Xr = 20.9 for Model CURV-RIM; see 
Tab. Eland Fig. [T^] and [12]) . We would like to point out that adding even larger dust grains 
cannot solve this discrepancy, since for grain sizes > 1.2 /im the inner r im location becom es 



practically independent of the grain size distribution (see discussion in llsella et al.ll2006l ) 



4.4. Models VERT-RIM-ACC & CURV-RIM-ACC: Dust disks with active 

inner gaseous accretion disk 

In passive circumstellar disks, the infrared emission is generally assumed to originate 
almost entirely from dust; the emissivity of the inner dust-free gaseous part of the disk, at 
radii smaller than the dust sublimation radius, is negligible. In an actively accreting disk, 
on the other hand, viscous dissipation of energy in the inner dust-free gaseous part of the 
accretion disk can heat the gas to high temperatures and give rise to significant amounts 
of infrared emission from optically thick gas. The inner edge of this gas accretion disk is 
expected to be located a few stellar radii above the stellar surface, where the hot gas is 
thought to be channeled towards the star via magnetospheric accretion columns. While the 
magnetospheric accretion columns are too small to be resolved in our interferometric data 
(3 -R^ correspond to 0.09 AU or 0.12 mas), the infrared emission from hot gas inside the 
dust sublimation radius should be clearly distinguishable from the thermal emission from 
the dusty disk due to the different temperatures of these components and the resulting 
characteristic slope in the NIR- and MIR-visibilities. 



As MWC 147 is a quite strong accretor (Mace ~ 10"^ M^v r~^ iHillenbrand et al.lll992r) 



infrared emission from the inner gaseous accretion disk is likely. iMuzeroUe et al.l (120041 ) found 
that even for smaller accretion rates the gaseous inner accretion disk is several times thinner 
than the puffed-up inner dust disk wall and is optically thick (both in radial as well as in 
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the vertical direction) 



To include the thermal emission from the inne r gaseo us disk in our radiative transfer 
models, we use a similar approach as lAkeson et al. J 20051). We assume that the accretion 



luminosity is emitted from a viscous accretion disk (jPringlelll98ll ) which emits at each radius 
r as a black-body of temperature 

4 , , (?,GM^m\ ( , -\i/2 

^sas(r) = (^^^ j (l - v/iV^) . (6) 

We run r from the magnetic truncation radius Rcomt to the dust sublimation radius. To 
estimate -Rnnrot for MW C 147, we use the measured rotation velocity {vsini = 90 kms^^. 



Boehm fc Catalalll995l ). the stellar parameters from Tab. [H and an intermediate inclination 
angle (as derived from the geometric model fits in Tab. El i = 60°), and yield -Rcorot = 
GM^/{27rv^) ~ 3R^. In our Monte Carlo radiative transfer simulation, the photons from the 
gaseous accretion disk are emitted isotropically from the disk midplane and then propagate 
through the simulation grid. In the course of a simulation, the accretion disk is assumed 
to be totally optically thick. Another possible consequence of the presence of an optically 
thick inner gas disk may be shielding of the inner dust rim, allowin g dust to exist c l oser t o 



the star. However, based on the results of the theoretical work by iMuzeroUe et al.l (l2004j ). 
we consider this to be a secondary effect and do not include it in our modeling. In order to 
reproduce the SED in the presence of the excess infrared continuum emission from optically 
thick gas, we vary the fraction between small and large dust grains in the disk in radial 
direction, which is in qualitative agreement with the indications for grain growth found in 
the inner disk regions (Sect. 



Including the accretion luminosity from an inner gaseous disk improves the agreement 
between model predictions and observed visibilities considerably. With a flared disk geometry 
and an accretion rate of Mace = 7 x 10~^ MQyr^^ (see Tab. E]), both the SED and the 
interferometric visibilities (Fig. [H] and [TH) are reasonably well reproduced. This result is 
not very sensitive to the precise shape of the inner rim (Model VERT-RIM- ACC Xr = 0.99; 
Model CURV-RIM-ACC: Xr = 1-24). 



5. Conclusions 

We have presented infrared long-baseline interferometric observations of MWC 147, 
constraining, for the first time, the inner circumstellar environment around a Herbig Be star 
over the wavelength range from 2 to 13 /im. 



-38- 



The interferometric data obtained from the PTI archive and with VLTI/ AMBER sug- 
gest a characteristic diameter of just ~ 1.3 AU (Gaussian FWHM) for the NIR emitting 
region, while the MIR structure is about a factor of 7 more extended (9 AU at 11 fim). 
Within the i^-band, we measure a significant increase of size with wavelength. Comparing 
the measured wavelength-dependence of the characteristic size with the class of analytic ac- 
cretion dis^jnod^lSj_wiin3h_are currently most commonly applied fo r the modeling of HAeBe 
SEDs (JHillenbrand et al.lll992l ) and NIR interferometric data (e.g. IVinkovic fc Jurkig 120071 ) 
yields that these models cannot well reproduce the measured significant increase of the appar- 
ent size with wavelength. To test whether more realistic physical models of the circumstellar 
dust environment yield better agreement, we employed 2-D radiative transfer modeling. The 
radiative transfer models were constructed to fit the SED from 0.3 to 450 /xm, including an 
extended envelope and an inner dust shell/disk geometry. While models of passive irradi- 
ated disks with or without puffed- up rims are able to reproduce the SED, they are in conflict 
with the interferometric observables, signiflcantly overestimating the size of both the NIR 
and MIR emission. This is also the case, to an even bigger extend, for spherical shell ge- 
ometries. Therefore, in our radiative transfer models, we incorporated additional accretion 
luminosity emitted from an inner gaseous disk, yielding significantly better agreement with 
the interferometric data. The best-fit was obtained with a fiared Keplerian disk seen under 
an inclination of ~ 50°, extending out to 100 AU and exhibiting a mass accretion rate of 
7 X 10"*^ MqYt~\ 

Since MWC 147 belongs to the group of "undersized" Herbig Be star disks (which 
have measured NIR diameters smaller than expected from the size-luminosity relation; see 
Monnier et al.ll2005l ). our detailed spectro-interferometric study confirms earlier speculations 
that the NIR emission of Herbig Be stars in this group might be dominated by gas emission 
from inside the dust sublimation radius. 



Furthermore, our study demonstrates that the spectro-interferometric capabilities of the 
latest generation of long-baseline interferometric instruments are particularly well suited to 
reveal the contributions from active accretion processes taking place close to the star. To 
achieve further progress on MWC 147, it seems promising for future spectro-interferometric 
observations not only to extend the spectral coverage (e.g. to the J- and if-band) but to 
employ also a higher spectral r esolution in order to stud y the spatial origin of accretion- 
related emission lines (e.g. Br7, iGarcia Lopez etal.l 120061), maybe trac ing ionized, optically 
thin gas in the disk or magnetospheric accretion (JHartmann et al.lll994l ). In addition, future 
long-baseline interferometric observations will measure the closure phase relation, which is a 
sensitive measure to asymmetries in the source brightness distribution, and which provides 
also the key to direct aperture synthesis imaging of the sub-AU circumstellar environment. 
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